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FINITE APPROXIMATIONS OF PHYSICAL MODELS OVER 

LOCAL FIELDS 


ERIK M. BAKKEN AND TROND DIGERNES 
Dedicated to Raja - dear friend, colleague, and inspirator. 


Abstract. We show that the Schrodinger operator associated with a physical 
system over a local field can be approximated in a very strong sense by finite 
Schrodinger operators. Some striking numerical results are included at the end 
of the article. 
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1. Introduction 

In |DVV94) it was shown that a quantum mechanical Hamiltonian of the form 
H = —A + V, acting in L^(R‘^), with potential V(a:) —>■ oo as |a:| —>■ oo, can be 
approximated in a very strong sense by finite quantum systems. In this note we 
present a similar theorem for quantum systems over a local field K. 

The results of [DVV94) were later extended to a setting of locally compact abelian 
groups in |AGK00j . The results of the latter thus supersede both those of [DVV94) 
and of this article. However, the proofs of |ACK00] used non-standard analysis. We 
have found it worthwhile to present a proof which does not rely on non-standard 
methods. 

In [DVV94) two proofs of the main convergence theorem were given: a func¬ 
tional analytic one and a probabilistic one. The latter gave a somewhat stronger 
convergence result for stochastic Hamiltonians. In the present note only functional 
analytic methods will be considered. A stochastic proof will be discussed in a 
forthcoming paper. 


The research of the second named author was partially supported by the Norwegian Research 
Council. 
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In an earlier article [BDLWl^ finite approximations over Qp were treated. The 
current article supersedes that one; also, the proofs which were omitted there, are 
given here. 

In Section [2] we give a quick review of local fields. In Section [3] we construct 
finite models for the Schrodinger operator over a local field, and in Section |T] we 
prove the main convergence theorem. 

In Section [S] we use our finite models to carry out a numerical investigation of 
the Schrodinger operator over the quadratic extension Q 3 [-\/ 3 ] of Q 3 . We show that 
there is remarkable agreement between numerical and theoretical values for both 
eigenvalues and eigenfunctions. Both types of eigenfunctions (radial ones and those 
supported on single shells) appear already at the finite level. 

2. Local fields 

We give here some quick facts about local fields. For a thorough treatment, see 
the classic treatise of A. Weil |Wei74[ Ch. I]; for a quicker review, see the book of 
Kochubei |KocOI[ Ch. 1.3]. 

A local field is a non-discrete, locally compact field. The only connected local 
fields are R and C. Disconnected local fields are, in fact, totally disconnected. 

Every local field comes equipped with a canonical absolute value which defines 
its topology. It is is induced by the Haar measure and is called module in [Wei74] . 
It is Archimedean in the case of R and C, and non-Archimedean in all other cases; 
it coincides with the usual absolute values for the fields R, C, and Qp. For a 
general local field K we will denote the canonical absolute value by | • | (or by | • jx 
if needed for clarity); for Qp we will denote it by | • |p. 

Convention. Since all local fields except R and C are (totally) disconnected, it is 
customary to reserve the term ’local field’ for a (totally) disconnected, non-discrete, 
locally compact field. We will follow that convention here. 

With this convention, there are two main types of local fields: 

Characteristic zero. The basic example of a local field of characteristic zero is the 
p-adic field Qp {p a prime number). Every local field of characterisitic zero is a 
finite extension of Qp for some p. 

Positive characteristic. Every local field of positive characteristic p is isomorphic 
to the field Fq((t)) of Laurent series over a finite field F,, where q = p^ for some 
positive integer / > 1 . 

Let AT be a local field with canonical absolute value | • |. Following standard 
notation, we set 

O = {cc G a: : |a;| < I}, P = {a; G AT : |a;| < 1}, U = 0\P. 

O is a compact subring of AT, called the ring of integers. It is a discrete valuation 
ring, i.e., a principal ideal domain with a unique maximal ideal. P is the unique 
non-zero maximal ideal of O, called the prime ideal, and any element /3 G P such 
that P = PO is called a uniformizer (or a prime element) of AT. For Qp one can 
choose P = p, and for Fq{{t)) one can take P = t. 

The set U coincides with the group of units of O. The quotient ring O/P is a finite 
field. If q = p^ is the number of elements in O/P {p: a prime number, /: a natural 
number) and /3 is a uniformizer, then |,5| = 1/q, and the range of values of | • | is 
q^, N G Z. Further, if S' is a complete set of representatives for the residue classes 
in O/P, every non-zero element x G K can be written uniquely in the form: 

X = P~'^{xq + XiP + X2P‘^ H-), 

where m G Z, Xj G S, xq ^ P. With x written in this form, we have jxj = 9 "*. 
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For a general field extension K/F we use the following standard notation: / = 
index of inertia, and e = ramification index. These are connected through the 
formula [K : F] = ef. If e = 1, the extension is unramified, and if / = 1, the 
extension is totally ramified. 

2.1. Characters and Fourier transform. We first fix a Haar measure fj, on K, 
normalized such that /r(0) = 1. The Fourier transform on if is given by 

(•^/)(C) = / f{x)x{-xOdx, 

JK 

where y is a suitably chosen non-trivial character on K, and dx refers to the Haar 
measure lust introduced. For our set-up it will be essential to use a character of 
rank zertfl. We describe a procedure for achieving this in the two main cases: 

2.1.1. Case 1: char K — 0. In this case if is a finite extension of Qp, and a character 
of rank zero is obtained by setting 

Xix) = Xp (Tric/Qp(/3"‘'a;)) , x G K, 


where 

• Xp is the canonical character on Qp—i.e., Xp(a:) = exp(27rf{a;}), {a;} = 
fractional part of x. 

• Tr^/q : if —Qp is the trace function associated with the extension 

K/ClJ. 

• P is a. uniformizer as defined above. 

• d is the exponent of the different of the extension K/C^p. It is the largest 
integer d such that Tr^/Q^(a:) G Zp for all x with |a;| < q‘^ (note that d > 0 
since Tr^^/q^ : O —^ Zp). 

2.1.2. Case 2: char if > 0. In this case we may identify if with the field Fq((t)) 
of Laurent series in the indeterminate t with coefficients from the hnite field F^, 
q = pd ^ consisting of elements of the form x = XiC, Xi G F^, m G Z. Let rj 

denote the canonical character on Fg, i.e., ri{x) = exp TrF_j/Fp(a;)j, and define 

X(a:) = ry(x_i), 

where x_i refers to the expansion x = YlTLm xf^- Then x is a rank zero character 
onif = Fg((i)). 

Notice that any Fourier transform based on a rank zero character is an L^- 
isometry with respect to the normalized Haar measure defined above (since F\o = 
lo for any such Fourier transform F', here and elsewhere 1 denotes characteristic 
function). Thus = F* is given by 

{F-^ f){x) = {F* f){x) = [ f{y)x{xy)dy. 

JK 

Convention. For the rest of this article F will denote a Fourier transform based on 
a rank zero character on if. 


^We remind the reader that the rank of a character x is defined as the largest integer r such 
that xlfl.. = 1- 
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3. Finite approximations over a local eield 

Our object of study is a version of the Schrodinger operator, defined for Qp in 
the book of Vladimirov, Volovich, Zelenov [VVZ94] . and generalized to an arbitrary 
local field K by Kochubei in |Koc01) : 

H = D°‘ + V, 

regarded as an operator in L'^{K). Here a > 00 , Z? = where {Qf){x) = 

\x\f{x) is the position operator, and is the Fourier transform on L"^{K). V 
(the potential) is multiplication by a radial function: (Vf)(x) = v{x)f{x), v{x) = 
?p(|x|) for some function w defined on [ 0 , oo). We assume v to be non-negative and 
continuous and that v{x) —>■ oo as |a;| —>■ oo. 

Due to a conflict of notation later in this article, we will use the symbol P for 
the differentiation operator (instead of D). With this notation we have 

H = P“ -h V . 

The operator H has been thoroughly analyzed (see |VVZ94) for K = Qp and 
[KocOl] for general K)-. It is self-adjoint on the domain {/ G L‘^{K) ■ P°‘f + Vf€ 
L^[K)}, has discrete spectrum, and all eigenvalues have finite multiplicity. Our 
next task is to set up a finite model for this operator. 

3.1. Finite model. Keep the above notation, i.e.: K is a local field, q = p -^ is the 
number of elements in the finite field O/P, ,5 is a uniformizer, and S' is a complete set 
of representatives for O/P. For each integer n set = /3~^0 = ball of radius q”. 
Then Bn is an open, additive subgroup of K. For n > 0 we set Gn = Bn/B-n- 
Then Gn is a finite group with elements. Since the subgroup i?_„ will appear 
quite frequently, we will often denote it by Hn, to emphasize its role as a subgroup. 
So Hn = B-n = /5"0 = ball of radius q“", and G„ = H-n/Hn- Each element of 
Gn has a unique representative of the form -|-a_„ +i/3-'^+i + --- + a_i/3-i + 

Oo -I- ai/3 -I-h -I- at G S. We denote this set by and call 

it the canonical set of representatives for G„; we also give it the group structure 
coming from its natural identification with G„. 

Let again p denote the normalized Haar measure on K (cfr. 12.11) . Since Hn is an 
open subgroup of K, we obtain a Haar measure p,„ on G„ = by setting 

/4n(^ T Hn^ — M(^ T Hn') — pi^Hn') — Q , for X Hn G Gn- 

So each ’’point” x + Hn of G„ has mass g”", and the total mass of G„ is 

g2n . q-n ^ ^ 

With this choice of Haar measure on G„ the mapping which sends the charac¬ 
teristic function of the point x Hn in G„ to the characteristic function of the 
subset x + Hn of K, is an isometric imbedding of Lf[Gn) into L‘^{K). We regard 
operators on Lf{Gn) as operators on Lf{K) via this imbedding, by setting them 
equal to 0 on the orthogonal complement of the image of L^(G„) in Lf{K). 

We introduce the following subspaces of Lf{K), along with their orthogonal 
projections : 

• Crt = {/ € L^(iL)| supp(/) C Bn}. The corresponding orthogonal projec¬ 
tion is denoted by G„ and is given by: Gnf = 

• Sn = {f & L‘^{K)\f is locally constant of index < g“"}. The correspond¬ 
ing orthogonal projection is denoted by Sn and is given by: 

{Snf){x) = q'^ Jh„ f{x+y)dy = f{x + y)dy = a.ye{f,n,x), where 


^For a direct analog of the Laplacian one should set a = 2. However, as is customary in the 
non-Archimedean setting, one works with an arbitrary o > 0, since the qualitative behavior of 
the operator H does not change with o > 0. 












FINITE APPROXIMATIONS 


5 


we have introduced the notation ave(/, n, x) for the average value of / over 

X + Hn. 

• Vn = C„ n The corresponding orthogonal projection is denoted by -D„. 
Note that L^(Gn) is mapped onto via the isometric imbedding mentioned above. 
Thus L‘^{Gn) can be thought of as the set of functions on K which have support in 
Bn and which are invariant under translation by elements of Hn (= B-„). 


Lemma 3.1. The projections Cn and Sn commute, thus the projeetion D„ onto 
the subspace is given by: 


Dn = CnSn = 


Proof. 

(5'„C'„/)(a;) = g” / {Cnf)ix + y) dy = [ 1 b,,{x + y)f{x + y) dy 

J Hrx JH-n 

=\" / lB„(a;)/(a; + 2 /)dj/= / f{x + y)dy 

Jn-n, JH-n 

= (GnSnf)(x) 

where the equality (*) follows from ultrametricity, namely: x + y G Bn x G Bn 
when y G Hn = B-n- D 

We next show that the Fourier transform behaves nicely with respect to these 
subspaces. 


Proposition 3.1. We have: 


Cn — Bn^ Sn — tjn, and hcnCe P n — n, 

and the same relations hold with in place of P. 

As a consequence, the following commutation relations hold: 


PCn = S'nJ", PSn = GnP, 


PDn = DnP. 


Proof. Let f G Cn and take any h G Hn. Then 

{J^f){^ + h)= [ f{x)x{-x{f + h))dx= [ f{x)x{-xf)x{-xh)dx 

Jk J Bn 

= [ fix)x{-xf)dx = [ f{x)x{-xf)dx 

JBr, JK 


= i^m) 


since \xh\ < g” • = 1 and y has rank zero. This proves PCn C 5„. 

Next let f G Sn and assume iPf){£,) ^ 0, |^| = g™. We must show that m < n. 
For any h with |ft,| < we have 

f{x)x(.-xOdx = [ f{x + h)x{-x^)dx 

JK JK 

= [ f{x)x(.-{x-h)0dx = [ f{x)x{-xOx{hf)dx 
JK JK 


which, since iPf)i£,) 0, gives 

x{d(f) = 1 for all h G i3_„ . 

This means that y is identically equal to 1 on the ball f • i3_„ = B-n+m, and since 
X has rank zero, we must have —n + m < 0, i.e., m < n. This proves PSn C C„. 
Since obviously the same relations hold with in place of P, we have equalities 
everywhere, i.e., PCn = Sn, PSn = C„, and hence PVn = Bn. 



6 


ERIK M. BAKKEN AND TROND DIGERNES 


As for the commutation relations: The relations just proved - and the same 
ones with = T* instead oi T - imply that SnTCn = = TSn, 

Sn^*Cn = ^*Cn, CnT*Sn = J^*Sn- Taking adjoiuts and combining, we get 
TCn = SnJ^ and TSn = CnT. Multiplying TCn = by Cn on the left and 
multiplying iFSn = by Cn on the right gives CnJ^Cn = CnSn^ = J^SnCn, 

i.e., J-Dn = DnT. □ 


3.2. Fourier transform at the finite level. We need to establish a relation 
between the Fourier transforms on K and Gn- 

So let as before y be a rank zero character on K and let T be the associated 
Fourier transform. Like any additive character on a field, y gives rise to a symmetric 
bi-character A on AT by setting X{x,y) = xi^v)- It descends to a bi-character on 
Gn = Bn/B-n, since if x' = x+h, y' = y+k with /i, fc G i3-„, then A (a:-I-h, y+k) = 
X{{x + h){y X k)) = x{xy)xixk)xihy)xihk) = x{xy) = X{x,y) (the arguments in 
the last three factors of the product all have absolute value < 1). So we can define 
a bi-character Xn on Gn by setting Xn(\x\, [?/]) = x{xy). Since X is non-degenerate 
on K, so is Xn on G„. Indeed, if a; G S„, |a;| = g™, and Xn{[x\, [y]) = y(a:y) = 1 for 
all y G Bn, then y = 1 on the ball x ■ Bn = Bn+m, which implies Bn+m C Bq = O 
since y has rank 0. But this means that m -|- n < 0, i.e., m < —n, and so a; G B-n, 
i.e., a; = 0 as an element of Gn = 

Setting Xn,[y]{[x\) = A’„([a:], [y]), it follows that the characters Xn,[y] exhaust all of 
Gn as [y] runs through Gn, i.e., the bi-character Xn implements the self-duality 
of the finite abelian group Gn- The canonical choice for an L^-isometric Fourier 
transform on Gn is then given by (recall that Gn has elements): 


iXnf){[x]) = ^= ^ fi[y])Xni-[x],[y]) 

= g"” fily])x{-xy), [x] GGn, / G L^{Gn) , 
[?/] 6 G„ 


(3.1) 


or, in terms of the set of representatives Xn, 


{Bnf){x) = q " ^ f{y)Xn{-x,y) 
vex„ 

= q~'^ f{y)x{-xy), xeXn, fGL‘^{Xn)- 

yex„ 


(3.2) 


The following result is now more or less obvious, but we state it as a proposition 
because of its importance. It plays a crucial role in the proof of the main convergence 
theorems, and simplifies matters considerably compared to the situation over R, 
where the relation between the finite and infinite Fourier transform was much more 
complicated (see |DVV94[ p. 626-627]). 

Proposition 3.2. Let the Fourier transforms F and Fn he as above. Then F 
leaves the space Vn (— L^(G„)) invariant, and 


F\x)n — ^‘I.e., Fn — FDn — L)nF . 


( 3 . 3 ) 
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Proof. The first part of the statement has already been proved (Proposition 13.II) . 
For the second part, take any f G Vn and let x G Bn. Then: 

{Bf){x)=[ f{y)x(.-xy)dy^ = " [ f{y)x{-xy)dy (3.4) 

JK J Br, 

= / fiy)x{-xy)dy= f{z)x{-xz)q-^ (3.5) 

dz+H„ 

= {Bnf){x), (3.6) 

where the equality {*) follows from the fact that the function y —)> f{y)x{~xy) is 
constant on z + Hn (since x G Bn) and y{z + Hn) = □ 


3.3. Dynamical operators at the finite level. For the finite versions of the 
dynamic operators we take their compressions by Dn, i.e., 14, = DnVDn, Qn = 
DnQDn, Pn = = DnT~^QBDn = T~^QnT = J^n^Qnd^n- Before Com¬ 

puting what these operators do to an / G L^(G„), let us find out what the projection 
Sn does to a radial function v(x) = r(;(|a:|): 


{Snv){x) 


1 

y{Hn) 



v{x + h) dh 


1 

y{Hn) 



w(\x -\- h\) dh 


ave(v, n, x) 


(*) 


v(x), 

ave(v, n, 0), 


|a:| > g-" 

\x\<q-^, 


where again ave{v,n,x) means the average value of v over x + Plm and where 
ultrametricity was used in the equality (*). 

Next we compute the effect of the finite operators on an / G L^(G„). For 14 we 
get, remembering that V is multiplication by a radial function v: 


{Vnf){x) = {DnVDnf){x) = {CnSnV f){x) = 1b„ (x) (S'n fo/) (cc) 

= iRn(a^) ,]r X / {Vf){x + h)dh 
h{Hn) JHn 

= 1 b„(x) / v{x + h)f{x + h)dh 

h\Hn) JHn 


= '^bAx) J / v{x + h)dh]f{x) 

h\Hn) JH„ 


= ave(w, n, x)f{x) 


v{x)f{x), 

ave(w,n,0)/(0), 


|x| > g-” 
|x| < g-” 


{Vf){x), \x\> q ” 

ave(w,n,0)/(0), |x| < g"" 


(/ G B{Gn)) . 


In particular, for the operator Qn this gives, writing g(x) = |x|: 


{Qnf){x) = ave{q,n,x)f{x) 


|x|/(x), \x\> q ” 

ave(g,n,0)/(0), |x| < g"" 


{Qf ){x), kl > g ” 

ave(g,n,0)/(0), |x| < g"” 


(/ G B{Gn)) . 


For Pn we get 

Pnf = DnPDnf = GnSnPf = GnSnB~^QPf = GnP-^GnQTf = GnP-^QBf 

= CnPf if G B{Gn)) . 
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We now set Hn — + Vn, the Hamiltonian for the finite model, and aim to show 

that the analog of Theorem 4 in |DVV94] holds in the present setting. 

4. Convergence of the finite models 

Keep the notation and assumptions of the previous section. 

There are two main steps to proving the analog of Theorem 4 of [DVV94| : Estab¬ 
lishing the convergence Hn —>■ H in the strong resolvent sense, and proving a form 
of uniform compactness for the resolvents (J -I- The proofs follow a pattern 

similar to that of [DVV94) . but we are able to simplify some of the arguments, 
partly due to the non-Archimedean nature of K. 

As for strong convergence of the resolvents: According to [KocOlj . Section 3.2, 
the space T) of locally constant functions with compact support is a core for the 
Hamiltonian H = P“ -|- V. Hence it is a common core for all the Hamiltonians 
Hn > 1) and H. For f £ V we have / G Vn for large n, hence lim^Q"/ = 
lim„ DnQ°‘f = Q°‘f in the strong operator topology; further: P“/ = T~^QnTnf = 
J^-^QnJ^f ^ = P“/,0 and Vnf ^ Vf. Thus = (P“ + H„)/ ^ 

(P“ -I- V)f = Hf for all f £V. Here we have used the obvious fact that Cn —>■ I, 
and hence Sn = TCn^~^ —>■ I and = CnSn = SnCn —>■ I, in the strong 
operator topology. 

By Theorem VHI.25 of |R,S80] it now follows that Hn —>■ H in the strong resolvent 
sense. 

The compactness of the resolvent (/ + H)~^ follows by classical arguments (see, 
e.g., [DVV94[ p. 623] for the case the same proof works for L?{K)). For the 

resolvents (/ -f we need a form of uniform compactness which is formulated 

as follows: 

Definition 1 (Uniform compactness). A sequence of bounded operators (M„) on a 
Hilbert space H is said to satisfy a condition of uniform compactness if the following 
conditions hold: 

(1) The sequence (Mn) is uniformly bounded. 

(2) There are subspaces with invariant under M„ such that for every 
sequence ((?„) with gn £ Ln and ||(7n|| < 1, the sequence {MnQn) is relatively 
compact in H. 

Remark. Notice that the individual operators are not required to be compact 
on H (and in our applications they will not be). Still, if the above conditions 
are fulfilled, we will say that the sequence (M„) is uniformly compact, even if the 
individual M„ are not compact. 

For our purposes the usefulness of uniform compactness lies in the following two 
results. They give a strong connection between the spectral data of the operators 
in an approximating sequence (M„) and their strong limit M . 

Lemma 4.1. Let M„, be as in Definitions^ and assume that the sequence 

Mn converges strongly to a bounded operator M. Assume further that there are 
eigenvectors gn and corresponding eigenvalues A„ such that gn £ Ln, Hffnll = 1 
and Mngn = ^nSn- Then any non-zero cluster point Aq of the sequence (Xn) is an 
eigenvalue of M, and there is a subsequence of (gn) which converges to a vector g 
such that Mg = X^g. 


^Proving the limit P^f —t P°‘ f required considerable effort in IDVV941 . due to the fact that 
the Fourier transforms at the finite and infinite level did not match up nicely. Here the finite 
Fourier transform is simply the restriction of the infinite one, and the limit becomes a triviality. 
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Proof. By uniform boundedness, all the A„ are confined to a bounded set. Hence 
there is a subsequence of (A^,) (still written (A„) after re-indexing) which converges 
to a scalar Aq, say, with Aq ^ 0. By uniform compactness, M„gn has a convergent 
subsequence (again written Mngn after re-indexing). It follows that the sequence 
gn = -^Mngn converges to an element g, say. Since —>• M strongly, it follows 

that Mngn Mg; indeed, remembering that the Mn are uniformly bounded: 
II Mngn - Mg II<11 Mn || || “ 5 II + II “ Mg ||-)> 0. So altogether we have: 

g = \imgn = lim = -^Mg, i.e., Mg = Xog. □ 

Notation: We let o'p{A) denote the set of positive eigenvalues of an operator A. 
Further, for a self-adjoint A we let denote the projection valued measure of A, 
and for a projection E we let r{E) denote its range. 

Proposition 4.1 (Cfr. Lemma 3 in [DVV94) 1. Keep the notation and assumptions 
of the previous lemma. In addition, assume the following: (i) The operators Mn, M 
are self-adjoint, and 0 < M, M„ < I, (ii) M is compact on TL, and Mn is compact 
on Ln. Then the following hold: 

(1) If J is a compact subset of (0,1] with J fl (jp{M) = 0, then J fl crp(M„) = 0 
for large n. 

(2) If X G ap{M), there exists a sequence (A„) with A„ G a{Mn) such that A„ —>■ 
A. Further, if J is a compact neighborhood of an eigenvalue X G ap{M), 
not containing any other eigenvalues of M, then any sequence (A„) with 
Xn G <jp{Mn) n J converges to X. 

(3) Let X and J he as in (2). Then dimP^"(J) = dimP^(J) for large n, and 
for each orthonormal basis {ei,..., Cm} for r [P^{j)) there is, for each n, 
an orthonormal basis {e",..., e((j} for r (P^"(J)) such that lim„_).oo e" = 
Ci, i = 1,... ,m. 

Proof. (1) If J n ap{Mn) ^ 0 for arbitrarily large n, there are infinitely many A„ 
in J. The sequence (A„) thus has a cluster point in J, and hence, by the previous 
lemma, M has an eigenvalue in J. 

(2) The first part follows from the fact that Mn —>■ M strongly |RS80l Thm. 
VIII.24, Vol. 1]. Now let (A„^) be all the eigenvalues of the various M„ which lie in 
J, indexed in an arbitrary fashion. Then {Xn^) has a cluster point in J, which by 
the previous lemma is an eigenvalue of M. Since M has exactly one eigenvalue in 
J, it follows that the sequence (An,j,) has exactly one cluster point in J, i.e., (A„;,) 
converges to A. 

(3) For ease of notation set En = P^^{J) and E = P^^{J). We first prove 
that dimi?„ < dimi? for large n. Assume otherwise, and set m = dimi?. Then 
there exists a subsequence En^. of En such that dimA^j, > m for all k. For each 
k, choose m -I- 1 orthonormal eigenvectors ,..., e!jn+i for r(i?„j,). By uniform 
compactness there is a subsequence of (e^) which converges to an eigenvector for 
M. Repeating the process for each of the remaining eigenvectors, we obtain a set of 
TO-|-1 orthonormal eigenvectors for M, a contradiction. This proves dim En < dim E 
for large n. The converse inequality dim En > dim E follows from |II,S80l Thm. 
VIII.24, Vol. 1]: Since M„ —>■ M strongly, then En —>■ E strongly. For finite 
dimensional projections this implies dimifn > dimi? for large n. 

For the last statement take any orthonormal basis {ei,...,em} for r{E). Let 
us first show that the set {i?„ei,..., is linearly independent for large n. 

Assume to the contrary that it is linearly dependent for arbitrarily large n, and let 
1 > e > 0 be given. By strong convergence there is an no such that || EnCj — Cj ||< e 
for n > no, j = 1.. .m. Pick an n > no such that the set ..., A„em} is 

linearly dependent. From a linear dependence relation for this set, pick the term 
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with the largest coefficient 


EnCi, say - and solve for it. Then we have 


m 

Eji^i — ^ ^ OijE^Gj 

i=ij#* 


with \aj \ < 1. Take the inner product with on both sides to get 

m 

{En^^iT ^ ^ O^ji^EyiGj^Bi) , 

i=lj/* 


(4.1) 


which gives 

m 

{En^i Ci, Ci) T 1 — ^ (^■^) 

i=lj/i 

For the left hand side we have \{Enei — ei, ei) + l| > 1 —e, and for the right hand side: 

I j=jti — ej,ei)| < (m — l)e. For e < Ifm this gives a contradiction. 

Hence the set ... ,Enem} is linearly independent for large n. Now perform 

a Gram-Schmidt orthonormalization on this set to obtain an orthonormal basis 
{e”,..., for r{En). An elementary, but somewhat tedious, calculation then 
shows that lim„_>.oo e" = Cj, i = 1,.. ., m. □ 


We are now ready to prove a key result, namely that the sequence (/ + 
is uniformly compact in the sense of Definition [T] (see Proposition 14.HI) . This will 
pave the way for establishing our main result iTheorem 14.111 . To prove uniform 
compactness we will use the following version of the Kolmogorov-Riesz compactness 
criterion; it is proved for the case L^(R‘^) in |HOH10l Corollary 7], and the same 
proof works for L^{K): 

Proposition 4.2. Let F be a subset of L^{K). Then F is relatively compact if the 
following conditions are fulfilled: 

(1) sup^g^ |1/||2< C for some positive constant C. 

(2) limr_>oo sup^g^ l/(a:)P dx = 0. 

(3) limp^oo sup/g^ /|^|>^ df, = 0. 

Proposition 4.3. With Mr, = {I + Hn)-^, Ln = Vr,~ L^{Gn), and n = L^{K), 
the resolvents (/ + E[n)~^ are uniformly compact in the sense of Definitions^ 


Proof. Let (g„) be as in Definition [T] and set /„ = (1 -f i7„) ^Qn- Then 

ll/n||2 +{HnfnJn) = ((/+ Hn)/n,/n) = (ffn,/n) < || II2 ||/n || 2 , 

1 /2 

and so || /„ || 2 < 1 since {Hnfn, fn) > 0 , and it also follows that || Vj /„ || 2 < 1 and 
II Pn^^fn || 2 < 1- The first of the last two inequalities gives: 


1 >||K)/Vn||^> 


|af|>r 


|/„(a;)F dx < 


' v{x)\fn{x)\'^ dx > inf v{x) 

\x\>r 

1 


|ai|>r 


\fn{x)\'^dx 


inf| 3 ,|>rw(x) 


0 . 


uniformly in n as r —>■ oo. For the inequality (*) we used that {Vnfn)ix) = v{x)fn{x) 
for |x| > r > q~^. 

Next we use the inequality || Pn^^fn || 2 < 1, valid for all n. First we note that 

II P^'^fn II 2 =11 l|2 = || ||2=|| Pfu lU 

= \\Sr,Q^/^Pfn\\2 
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which gives, for all p > 0: 

l>||^„Q“/'^/n|| 2 > f [ \iQ^/^Tfn)ix)\‘^dx 

J|a:|>p '^\x\'>p 

= [ |xr-|(^/n)(x)| 2 dx>p“ / \{Tfn)ix)fdx 

J\x\>o J\x\>o 


'\x\>p 


\{TJn){x)\^ dx < — -^ 0 


uniformly in n as p —>■ oo. For the equality (*) we used that is locally 

constant away from the origin. Uniform compactness of the (J + Hn)~^ now follows 

from Proposition 021 D 


It now follows that with Mn = (/ + Hn)~^ and M = {I+ H)~^ all the conditions 
of Proposition 14.11 are satisfied, and via spectral mapping we can state the analog 
of Proposition im for Hn and H: 

Theorem 4.1 (cfr. Theorem 4 in [DVV94) 1. (1) If J is a compact subset of 

[0, oo) with J n cr{H) = 0, then J fl a{Hn) = 0 for large n. 

(2) If X € o'{H), there exists a sequence (A„) with Xn G a{IIn) such that 
A„ —^ A. Further, if J is a compact neighborhood of an eigenvalue X G 
cr{H), not containing any other eigenvalues of H , then any sequence A„ 
with Xn G cr(i/„) n J converges to X. 

(3) Let X and J he as in (2). Then dimP'^"(J) = dimP^(J) for large n, and 
for each orthonormal basis {ei,..., Cm} for r (^P^{J)) there is, for each n, 
an orthonormal basis {e",..., elf^} for r ( J)) such that hm„_>.oo e" = 
Bi, i = 1,... ,m. 


5. Numerical investigation of the Schrodinger operator over Q3[-\/3] 

5.1. Overview. In |VVZ94l Ch. 3, Section XII] a detailed analysis was carried out 
on the spectrum of the p-adic Schrodinger operator, and in |Koc011 Ch. 3] a similar 
analysis was performed on the Schrodinger operator over a general local field. 

Let as before H = P“ + V denote the Schrodinger operator over a local field K. 
The eigenfunctions of H can be divided into two main types, corresponding to two 
complementary subspaces of L^{K): those which are supported on a single spherical 
shell (which we shall call shell functions), and those which are radiafl. Of these, only 
the shell functions are completely understood: They belong to eigenvalues which 
can be determined from Diophantine equations, and there are explicit formulae for 
them. For radial eigenfunctions no such explicit formulae seem to be known. 

In this numerical study we specialize to the case of the Schrodinger operator 
H = + Q^) of the harmonic oscillator over the local field Q3[-\/3], which is a 

quadratic and totally ramified extension of Q 3 . We were interested in the following 
questions: 

• Do eigenfunctions of both types (shell functions and radial functions) show 
up already at the finite level? 

• Is there good agreement between the theoretical and numerical eigenvalues? 

• Is there good agreement between the theoretical and numerical eigenfunc¬ 
tions? 

• Are multiplicities correct? 

The answer to all these questions was ’y^s’. To illustrate this, we sum up some of 
the results in Tabled! 

^With notation as in IVVZ941 IKocOll , the set of shell functions comprises all the type I func¬ 
tions plus the shell functions of type II; the radial functions are all of type 11. 
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5.2. More details about the numerical experiment. The extension Q3 [a/3 ]/Q3 
is totally ramified, so with notation as in section[ 2 ]we have e = 2, and hence f — 1 
since ef = [Q 3 [-\/ 3 ] : Qs] = 2. Further, from q = follows q = p = 3, and as 
uniformizer we can take (3 = \/3, hence \I3\ = 1/q — 1/3. For the exponent of 
the different we have d = 1, so the character y defined in subsection 12.11 becomes 

X{x) = exp (27ri{TrQ^[^]/p^(^ x e QsiVS]. 

For the finite model we did experiments with n = 1, 2, 3, 4, so we were working 
with finite grids of sizes |Xi| = 9, \X 2 \ = 9^ = 81, IX 3 I = 9^ = 729, and IX 4 I = 
9"^ = 6561, respectively. Of particular interest to us was how the eigenfunctions 
came out: Would they clearly exhibit characteristics as shell functions or radial 
functions? They did. To illustrate this we give in Table [5] an excerpt from the 
value tables of three eigenfunctions: one is radial, one is a linear combination of 
two shell functions, and one is a pure shell function. We also wanted to compare 
our numerically computed eigenfunctions to the theoretical ones (evaluated on the 
grid). To do this, we measured the distance from each of the former to the linear 
span of the latter. Up to machine accuracy (10“^®), the distance came out as zero. 
We find this quite remarkable. 
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Appendix A. Tables for numerical eigenvalues and eigenfunctions 

The tables in this section should be self-explanatorjU. The data are taken from 
a computer run with n = 2 (i.e., 81 points in the finite grid). Each of the fnnctions 
in Table [5] is represented with 28 values, with values coming from each of the 5 
shells which occur for n = 2. 


®In the estimate for the lowest eigenvalue in Tabled (first entry in column 1) we are assuming 
that the estimate given in IVVZ941 p. 190] is valid also in our setting. We havent checked this in 
detail, but there are strong indications that it is true. 


Table 1 . Numerical approximations to the spectral data oi H = 
over Qaiv^]- 


Theoretical 

eigenvalue 

Numerical 

eigenvalue 

Theoretical 

multiplic¬ 

ity 

Numerical 

multiplic¬ 

ity 

Type of 
eigenfunc¬ 
tion 

Comment 

0 < Ao < 9/13 
Pi 0.6923 

0.6684 

1 

1 

radial 


7 

4.6922 

7 

1 

radial 


7 

4.7158 

7 

1 

radial 


5 

5.0000 

2 

2 

shell function 

2 = 1 + 1: 
Coming 
from two 
different 
shells. 

9 

9.0000 

4 

4 

shell function 

All 

supported 
on the 

same 

shell. 

7 

40.5213 

7 

2 

radial 


40 + 5/9 = 
40.5555... 

40.5555 

2 

2 

shell function 

2 = 1 + 1: 
Coming 
from two 
different 
shells. 

41 

41.0000 

8 

8 

shell function 

8 = 4 + 4: 
Coming 
from two 
different 
shells. 

45 

45.0000 

24 

24 

shell function 

24 = 

12 + 12: 
Coming 
from two 
different 
shells._ 
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Table 2. Eigenfunctions for three different eigenvalues, 28 values 
for each function, coming from all the 5 shells. Both kinds of 
eigenfunctions occur (shell functions and radial functions). - Shell 
no. fc (fc = 2 , 1 , 0 , — 1 , —oo) is the shell |a;| = 3^ (so shell no. —oo is 
the shell \x\ = = 0 ). 




Eigenfunction for A = 5 




bigenmnction tor the lowest 

Eigenfunctions here are 

lin- 



eigenvalue A ~ 0.6684 . 

ear combinations of shell 



It exhibits a perfect radial 
behavior. Notice also that 

functions from two differ¬ 
ent shells (shells 1 and 0 ). 

As should be expected, the 
function below exhibits non- 
radial behavior, being non¬ 
constant on each shell where 
it doesn’t vanish (shells 1 
and 0 ). 



Eigenfunction for A = 9 

. It 

the function is strictly posi¬ 
tive, in accordance with the 
corresponding statement for 
the case K = in |VVZ94l 

p. 186]. 

exhibits a perfect shell func¬ 
tion behavior, with support 
on shell no. 1 . 


Shell no. 


Shell no. 


Shell no. 

3.5818432-10"^ 

— CX) 

1.8757870-10-1'’ Ri 0 

— OO 

-3.8765003- lO-'® Ri 0 

— OO 

5.5430722-10"^ 

2 

2.0896995 - lO-i'’ Ri 0 

2 

1.6021680- 10 - 1 ® Ri 0 

2 

5.5430722- lO-'^ 

2 

8.7737711 - lO-ii' « 0 

2 

-9.1411700- 10-11' Ri 0 

2 

5.5430722-10"® 

2 

-1.4801152- 10-1® Ri 0 

2 

5.1268297- lO-H' Ri 0 

2 

5.5430722 -lO-^ 

2 

3.0773313-10-1® Ri 0 

2 

2.7677667-10-1® Ri 0 

2 

5.5430722 -lO-^ 

2 

-4.5409159- 10-11' Ri 0 

2 

-4.5822760- lO-i® Ri 0 

2 

5.5430722 -lO-^ 

2 

-1.0479409- 10-1® Ri 0 

2 

-1.3758518-10-1® Ri 0 

2 

5.5430722 -lO-^ 

2 

-2.3471948- lO-H" Ri 0 

2 

2.1385872- lO-H' Ri 0 

2 

5.5430722 -lO-^ 

2 

7.9466194- lO-H' Ri 0 

2 

-1.0549816- 10-1® Ri 0 

2 

5.5430722- lO-'^ 

2 

2.3950293- lO-i® Ri 0 

2 

2.3917324- lO-i® Ri 0 

2 

5.5430722-10-^ 

2 

6.4773691 - lO-H' Ri 0 

2 

1.2912546- 10-1® Ri 0 

2 

5.5430722- lO-'^ 

2 

-1.1431061- 10-1® Ri 0 

2 

-6.0210598- 10-11' Ri 0 

2 

5.5430722-10-^ 

2 

-1.3177515-10-11' Ri 0 

2 

-3.9251100- 10-11' Ri 0 

2 

5.5430722-10-^ 

2 

1.3595786-10-1® Ri 0 

2 

-5.0103544- 10-11' Ri 0 

2 

5.5430722 -lO-^ 

2 

3.2839452 - lO-H' Ri 0 

2 

1.2137971- 10-1® Ri 0 

2 

5.5430722 -lO-^ 

2 

7.8206625 - lO-H' Ri 0 

2 

-1.0063910-10-1® Ri 0 

2 

5.5430722 -lO-^ 

2 

3.3933100-10-11' Ri 0 

2 

-7.7900493-10-11' Ri 0 

2 

5.5430722 -lO-^ 

2 

8.8459742 - lO-H' Ri 0 

2 

2.2672330- lO-i® Ri 0 

2 

5.5430722 -lO-^ 

2 

2.2115193- 10-11' Ri 0 

2 

-1.1819127- 10-1® Ri 0 

2 

1.2747433-10"^ 

1 

-2.3459638-10-1 

1 

5.9907185-10-1' 

1 

1.2747433-lO-^ 

1 

2.3459638-10-1 

1 

-4.1084268- 10-1 

1 

1.2747433- lO-^ 

1 

-2.3459638-10-1 

1 

-1.0595734- 10-1 

1 

1.2747433-lO-^ 

1 

2.3459638-10-1 

1 

2.7644342 - 10-^ 

1 

1.2747433-lO-^ 

1 

-2.3459638-10-1 

1 

4.6050157-10-2 

1 

1.2747433-lO-^ 

1 

2.3459638 -lO-i 

1 

3.8319834-10-1 

1 

3.1960943-10-^ 

0 

3.9500330 -10-1^ 

0 

1.2637350-10-11 Ri 0 

0 

3.1960943-10-1 

0 

-3.9500330-10-2 

0 

-1.6035100-10-11' Ri 0 

0 

3.5768544-10-1 

-1 

2.2996138- lO-^ Ri 0 

-1 

-9.9411507- lO-ii Ri 0 

-1 


Department of Mathematical Sciences, The Norwegian University of Science and 
Technology, 7491 Trondheim, Norway 
E-mail address: erikmcLki@inath.iitnu.no 

Department of Mathematical Sciences, The Norwegian University of Science and 
Technology, 7491 Trondheim, Norway 
E-mail address: digernes@math.ntnu.no 
































